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Abstract 

We provide an integral formula for the free energy of the two-matrix model with polynomial 
potentials of arbitrary degree (or formal power series) . This is known to coincide with the r- 
function of the dispersionless two-dimensional Toda hierarchy. The formula generalizes the 
case of conformal maps of Jordan curves studied by Kostov, Krichever, Mineev-Weinstein, 
Wiegmann, Zabrodin and separately Takhtajan. Finally we generalize the formula found in 
genus zero to the case of spectral curves of arbitrary genus with certain fixed data. 



1 Introduction 

Many instances of integrable systems are obtained by means of a suitable limit (dispersionless or semi-classical) of 
a statistical theory. The departing point of our analysis in this paper is the random 2-matrix model [211 |U] , which 
is attracting growing attention due to its applications to solid state physics |15) (e.g., conduction in mesoscopic 
devices, quantum chaos and, lately, crystal growth[2*2*]'). particle physics |29| . 2<i-quantum gravity and string theory 
|HU13l l3*|. The model under inspection consists of two Hcrmitian matrices Mi, M 2 of size N x N with a probability 
distribution given by the formula 



d/i(M 1 ,M 2 ) = — &MidM 2 exp 
Zn 



1 Tr (Vi(Mi) + V 2 (M 2 ) - M X M 2 ) 



H 



K=l J=l 



where Vi are formal power series but soon will be restricted to polynomials for simplicity. The partition function 
Zn is known to be a r-function for the KP hierarchy in each set of deformation parameters (coefficients of Vi 
or V 2 ) and to provide solutions of the two-Toda hierarchy |28l ^ 12] < This model has been previously investi- 
gated in the series of paper ED E] where a duality of spectral curves and differential systems for the relevant 
biorthogonal polynomials has been unveiled and analyzed in the case of polynomial potentials. In [S] the mixed 
correlation functions of the model (traces of powers of the two non-commuting matrices) have been reduced to a 
formal Frcdholm-likc determinant without any assumption on the nature of the potentials and using the recursion 
coefficients for the biorthogonal polynomials. We briefly recall that the biorthogonal polynomials are two sequences 
of monic polynomials ([5] and references therein) 

ir n (x) = x n + ---, a n (y)=y n + ---, n = 0,l,... (1-2) 

that are "orthogonal" (better say "dual") w.r.t. to the coupled measure on the product space 

[dxdy TT n (x)a m {y)e~^ v ^ +v ^' x ^ = h n S mn , h n ^ Vn e N (1-3) 
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where V\ (x) and V2 (y) are the functions (called potentials) appearing in the two-matrix model measure . It is 
convenient to introduce the associated quasipolynomials defined by the formulas 

4> n {x) := - 7 ^ 7 r n _ 1 (x)e-*^W (1-4) 

My) ■■= -^vn-itoe-W* ■ (1-5) 

In terms of these two sequences of quasipolynomials the multiplications by x and y respectively are represented by 
semiinfinite square matrices Q = [QyjijeN* arL d P = [Pj]i jeN* according to the formulae 

Xlp n (x) = 2J Qn,mi>m(x) ; y(f>n(y) = P m , n (j) m {y) 
m m 

Q n ,m = = Pm,nj if n > m + 1. (1-6) 

The matrices P and Q have a rich structure and satisfy the "string equation" 

[P,Q] = hi . (1-7) 

We refer for further details to 01 El where these models are studied especially in the case of polynomial 
potentials. We also point out that the model can easily be generalized to accommodate contours of integration 
other than the real axes 0]|5]- 

The partition function is believed to have a large N expansion according to the formula 

--Lln^=^ = ^°) + ^W+... . (1-8) 

This expansion in powers of N~ 2 has been repeatedly advocated for the 2-matrix model on the basis of physical 
arguments |Kil 114) and has been rigorously proven in the one-matrix model |16j . In the two-matrix model this 
expansion is believed to generate 2-dimensional statistical models of surfaces triangulated with ribbon-graphs 
[131 111! IT7| . where the powers of iV -1 are the Euler characteristics of the surfaces being tessellated. From this 
point of view the term corresponds to a genus tessellation and the next to a genus one tessellation. 

The object of this paper is the leading term of the free energy, . It is the generating function of the 
expectations of the powers of the two matrices in the model 

(M\ K ) = Kd UK T^ + 0{N- 2 ) , (M 2 K ) = JdvjjM + 0{N- 2 ) . (1-9) 

Integral formulas for these partial derivatives at the leading order are known and involve integrals over a certain 
spectral curve (see below |18|V however a closed formula for the function J-^> itself was so far missing; this paper 
fills the gap (Thm. 12. 1|) . Remarkably, an algorithm for the computation of the subleading terms is known and also 
a closed expression of the genus 1 correction and therefore this paper precedes logically and complements 

|14j . The paper is organized as follows: in section^ we recall the main formulas known in the literature and set up 
the notation, linking our result with the relevant other approaches E3J EU ESI E3 E01 EDI CHI EH EE] The core of 
the paper is section |2~TI where the formula for the leading term of the free energy (dToda tau function) is presented 
and proved (Thm. I2.1fl . In section |2~21 we apply our result to the generating function of the canonical change of 
coordinates represented by the Lax operators of the dToda hierarchy. 

Finally in section |3 we extend the result obtained in section 12.11 and find an integral expression for the free 
energy of the two matrix model in the case the spectral curve is of arbitrary genus (Thm. I3.1fl . 



2 Planar limit (dToda hierarchy) 

In this section we investigate the planar free energy (J 7 ^ ' in the notation of the introduction) and we will make 
soon the common "one-cut" assumption (to be lifted in section|31) which amounts to saying that the multiplication 



2 



operators tend to meromorphic functions over a spectral curve of genus zero. Indeed in this limit the two multipli- 
cation operators for the wave-vectors defined by the biorthogonal quasipolynomials become commuting functions 
in the shift operator here replaced by the variable A ([7] and references therein) and they corresponds to the Lax 
operators of the dToda hierarchy: 

oo 

Q(A)= 7 A + X>A-i (2-1) 



fe=0 



3=0 

The original noncommutativity of the operators now translates to the following Poisson-bracket which is nothing 
but the dispersionless form of the string equation (|l-7|l 221 ED HH1 

fPO\-X dP9Q \ 9QdP -A (2 1\ 

{P,Q}.-A— — -A— — -1, (2-3) 

where t = HN in the large limit. The parameter t is the (scaled) number of eigenvalues of the matrix model and 
enters the relations 

±-j >P &Q = t- J-^QdP = t; (2-4) 

The contour chosen is a contour on the physical sheet of the Q (P respectively) plane around infinity (i.e. around 
A = oo, A = respectively). 

The deformation equations describe the infinitesimal variations of the operators P, Q under variations of the 
parameters of the potentials; they are known in the finite N regime as well (00 and references therein) whereas 
in the dispersionless limit are given by the evolution equations [231 EH ESI HZ] 

(d UK Q)x = {Q,(Q K )+o} ; (d Vj Q)x = {Q,(P J )-o} (2-5) 
{d UK P) x = {P, (Q K )+o} ; (pvj P)\ = {P, (P J )-o} (2-6) 

where the subscript o± denotes the positive (negative) part of the Laurent polynomial plus half the part constant 
in A, viz, e.g. 

(Q K U(\) = (Q K ) + (\) + \(Q K ) . (2-7) 
If the potentials Vt are polynomials of degrees di + 1, i = 1, 2 then both P and Q are finite Laurent polynomials 

Q(X)=~/\ + ^2a j \- j (2-8) 

fc=0 
di 

P(\) = l + J2&V ■ (2-9) 

In what follows we restrict to this case so as to avoid complication of convergence; however one may replace the 
contour integrals that will follow with formal residues of formal Laurent series and carry out the same computations. 
Another reason why we prefer the truncated setting is that then the two functions P, Q define a (singular) spectral 
curve of genus 5 = which is given by the polynomial locus (resultant-like) obtained from the determinant of the 
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following Sylvester matrix ([B], thanks to a remark by J. Hurtubise) 



= E(P,Q) 



d\+d 2 
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(2-10) 



This spectral curve is precisely the (limit of the) spectral curve of the four finite-dimensional folded differential 
systems for the quasipolynomials (SJ|7|. If we worked with formal power series it is not known whether a spectral 
curve in this sense can be defined. 

Taking this viewpoint A 6 CP 1 is the uniformizing parameter of the spectral curve E(P, Q). The two potentials 
are related to the parameters 7, ay, Pj by the relations 



P 

Q 



V{{Q) 



t 

p 



0(Q- 2 ) 

o{p- 2 ) 



near 00 p 



(2-11) 
(2-12) 



where the point ooq (oop) is the point on the spectral curve where Q (P respectively) has a simple pole; in the 
uniformization provided by the coordinate A it corresponds to A = 00 (A — respectively). By expanding both 
sides in powers of A and matching the coefficients one can realize that the coefficients of the two potentials are 
rational functions of the parameters 7, a 3 , Pj. More explicitly we have 



vi( q ) = y^ 



d ± + l 
K=l 



V 2 ( P ) 



d 2 + l 
,7=1 J 



U K 



1 

2iTT 



p 

7W 



dQ 



vj 



1 

2in 



Q_ 

P j 



dP 



(2-13) 
(2-14) 



The leading term of the free energy of the model is then defined by the differential equations 

dT „ 1 f „ „ K . „ 1 



du 



U K :-- 



K 



K2iir 



PQ^dQ 



T , res PQ ft dQ 



|^ = Vj := -J- / Q P J dP = - T res QP J dP . 

OVJ J ZlTT J J oo P 



(2-15) 



These equations are precisely the same that define the r-function of the dToda hierarchy. In the relevant 
literature |231 1241 |2*B1 127j the functions P, Q are the Lax operators denoted by L, L or C, C and the normalization 
is slightly different. 

We should also remark the following link to the works |2(J1 1191 13UI 1311 l2l)| in that if we take V\ = V = V2 we 
then have 7 £ R, aj — Pj. The function Q(A) is then the uniformizing map of a Jordan curve T in the Q-plane (at 
least for suitable ranges of the parameters) which is defined by either of the following relations 

P(1/X)=Q(X) ;|A| = 1 . (2-16) 

In the setting of |30U19ll2D| the function Q is denoted by z (and A by w) so that then P is nothing but the Schwartz 
function of the curve T, defined by 

z = s{z) , z e r . (2-17) 
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The coefficients of the potential V(x) = J2k T 2 ^ are ^ ne s °- caue d "exterior harmonic moments" of the region T> 
enclosed by the curve T 

t K = — I zz~ K dz (2-18) 

2«7T J r 

By writing z = S(z(w)) these integrals become residues in the w-plane. In our general situation the representation 
of the coefficients uk , v,j (|2-14|l cannot be translated to a surface integral, hence the need for a separate analysis. For 
conformal maps (i.e. Jordan curves) the r-function has been defined in |2(i| and given an appealing interpretation 
as (exponential of the Legendre transform of) the electrostatic potential of a uniform 2-dimensional distribution of 
charge in T> ^U] 



ln(rr) = - 







i 


1 






z 


~z~' 



d 2 zd 2 z' . (2-20) 



It can be rewritten as a (formal) series in the exterior and interior moments as 

21n(r r ) = -— ff d 2 z\z\ 2 + t w + ^ (t K w K + t K w K ) (2-21) 



47T_ K>Q 

where the interior moments arc defined by (the normalization here differs slightly from |3l)p 

w K = —[[ z K d 2 z , K>0 (2-22) 
tt-K JJt> 

w = ^ JJ ln|z| 2 d 2 z . (2-23) 

Unfortunately this formula is not exportable to our more general setting (in particular the logarithmic moment 
above) and to the general setting of the dispersionless Toda hierarchy; our principal objective is to fill this gap. 



2.1 Free energy of the 2-matrix model in the planar limit 

Let us focus on the function P as a (multivalued) function of Q (similar argument can be reversed for Q as a 
function of P); on the physical sheet P(Q) will have in general some branch-cuts with square-root singularities at 
the endpoint of each cut. These cuts are bounded in the physical sheet because P(Q) is analytic in a neighborhood 
of ooq. Note that 

V ^ = ^J^{ l -h) PdQ (2 " 24) 
V 2 (p) = ^Iln(l- |) QdP , (2-25) 

as one can immediately realize by expanding in powers of q, p. The integrals are well defined provided that q (p 
respectively) arc kept inside the contour of integration. In the following we will develop all the necessary arguments 
only for V\(q) and related objects, where the reader can obtain the relevant proofs for Vzip) by interchanging the 
role of P and Q. 

Let us now introduce the exterior potentials 2 

^i(q ut) = -^- —) PdQ (2-26) 

2 It has been pointed out to me by B. Eynard that they are sometimes referred to as effective potentials because they can be thought 
of as the effective potential felt by one eigenvalue in the field of the others. 
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In this formula the contour of integration is such as to leave the point q ou t in the outside region (whence the 
subscript). In general we will distinguish the choice of the contours by a subscript q out or qi n in what follows. 
The coefficients of $>i(q ou t) in inverse powers of q ou t are precisely (minus) the Uk defined in i|2-15|) . Note that 
expanding eq. (|2-26|) in inverse powers of q ou t the first d\ + 1 coefficients are the Uk coefficients as defined in 
(l2~T5l) . 

The first objective is to analytically continue $1 to the physical sheet so as to obtain a different representation 
of it. For this purpose we compute 

mout) = ±-i —7^ rPdQ = 1^- + -L / - 1 PdQ = (2-27) 

2«7r J q ou t (Q ~ qout) qout^in 2in J (Q - q out ) 

' ' ' ' -PdQ + P(q)= t ~-V((q)+P(q) , (2-28) 



q 2in J Q - q in " q 

where we have dropped the subscript outside of the integral as those terms are analytic functions in the whole 
physical sheet. Integrating once we obtain 

$i(Q) = -Vi(Q) + t ln(Q) + [ PdQ , (2-29) 

Jx q 

where X q is a point defined implicitly by the requirement <f>i = 0(Q^ 1 ) near ooq(= (A = oo)). Note that, in spite 
of the ln(Q) term, this is an analytic function around ooq. 
By similar reasoning we get 

$2(P) = ~ i / ln ( X ~ PJ ® dP = _ ^ 2(P) + 1 ln(P) + J x QdP (2_30) 

We need to introduce two more points (beside X q and X p ) on the spectral curve and a lemma: those are the points 
Aq, A p defined implicitly by the relations 

/ PdQ = Fi(Q(A)) >0 + t ln(A) + OiX- 1 ) , near oo Q (2-31) 

/ QdP = V 2 (P(A)) <0 + t ln(A) + 0(A) , near oo P . (2-32) 

By a simple inspection of the A°-coefficient in the LHS and RHS one finds that 

/ PdQ = f PdQ + [ PdQ = [ PdQ + 04(Q)) - t ln( 7 ) (2-33) 

JX q JA q JX q JA q 

[ QdP= [ QdP+ [ QdP= f QdP + (V 2 (P)) - t ln{j) , (2-34) 

Jx p JA P Jx p J A p 

where the subscript o denotes the constant part in A (which can be written as a residue). 
We now have the 

Lemma 2.1 The following relation holds 

li := Q(X P )P(X P ) + ( X PAQ = Q(X q )P(X q ) + [ X QdP = (pQ - Vi(Q) - V 2 (P)) + 2t ln( 7 ) . (2-35) 

Remark 2.1 The quantity fi will be proved in Corollarv l2.2l to be the derivative of the free energy w.r.t t, i.e. the chemical 
potential. It therefore corresponds to the logarithmic moment in the setting of |30l l2Q|[TIfll31| . 



G 



Proof. The equivalence of the two integrals is immediate by an integration by parts. And integration by parts is 
in fact the key to prove also the last part: 

J PdQ =PQ- P(A q )Q(A q ) - I QdP = PQ- P(A q )Q(A q ) - f QdP - f QdP . (2-36) 

J Ag J Ag JA q JA p 

Looking at the constant part in A in both sides of the above equation we conclude that 

(QP)o = P(A q )Q(A q ) + [ QdP . (2-37) 

Therefore we have 

(QP) = P(A q )Q(A q ) + [ QdP = P{A q )Q(A q ) + [ QdP + [ QdP = (2-38) 

JA q JX P JA q 

A X 

= P(A q )Q(A q )+ [ QdP + Q(X P )P(X P ) - Q(A q )P{A q ) - / PdQ = (2-39) 

JX V JA q 

QdP + Q(X P )P{X P ) - f PdQ - f PdQ = (2-40) 

'p J X q J Ag 

(Vi(Q) + V 2 (P)) - 2t ln( 7 ) + Q(X p )P(X p ) + f PdQ . (2-41) 



This concludes the proof of the Lemma. Q.E.D. 

We are now in a position to formulate the first main result of this paper. 

Theorem 2.1 The Free energy is given by the formula (up to constant) 

i — ZDV^ , . ™ {Vi(Q) + V 2 (P)-PQ, 



2T = res (P$idQ)+res (Q$ 2 dP)+- res (P 2 QdQ) + t res ^ v ; ^-dA +2t 2 ln( 7 ) (2-42) 

OOQ OOp 2 ooq v ' A=oo V A / 

= V u K U K + V vjVj + i res P 2 QdQ + t res ( Vi(Q) + V 2 (P) - PQ \ + ^ 



2 ooq A— oo \_ A 

Before proceeding to the proof a corollary and some remarks are in order. First off from the expressions in Thm. 
12.11 we find the well-known scaling property of the free energy ^3 1311 ITU| 

Corollary 2.1 The free energy defined in Thm. 12. II satisfies the scaling equation 

4^ = -t 2 + \ 2td t + ]T(2 - K)u K d UK + ^(2 - J)vjd Vj )T . (2-44) 

(More general scaling equations will be introduced later in Corollary 13. 2fl . 

Proof. The proof in the context of the normal matrix model can be found in the references quoted above and, in 
view of the formal equivalence of the normal matrix model with the two-matrix model |19j , the statement would 
follow also in our case. 

In order to be self contained we rederive this property in the present context; one way of proving formula 1)2-44(1 is 
from the expression (|2-43l) for T given in Thm. 12.11 by computing the residues involved after rewriting symmetrically 

the term res P 2 QdQ as i ( res P 2 QdQ + res Q 2 PdP I and computing it. A second, possibly instructive way to 

OOQ Z yOOQ OOp J 

derive it also directly from eq. (|2-42|l is by using the scaling properties of the various quantities involved. To this 
purpose we introduce the rescaling according to the formulas 

Q = SQ , P = 5P . (2-45) 



7 



Under this change of the functions Q,P we have (using formulas 12-41 12-15)1 1 

u K = S 2 ~ K u K ; vj = 6 2 ~ J vj ; t = 5 2 t 
Moreover, computing explicitly the residue res PQ'dX we find that 



A— oo 



7 2 = -t + Y^jocjPj 



(2-46) 



(2-47) 



from which one immediately obtains that 7 = $7. The exterior potentials are also conformally invariant, indeed 
$1 = -V X (Q) +t\n(Q) + I PdQ = ~5 2 Vi{Q) + 5 2 t\n{Q) + f PdQ + S 2 tln(5) + f PdQ (2-48) 

J X a J X a J X„ 



S 2 ®! + 8 2 tln{5) + [ PdQ 



(2-49) 



The last two terms cancel each other out because both the LHS and RHS must be 0(Q 1 ) = 0(Q 1 ). Repeating 
the argument for $2 we finally get 

$1 = <5 2 $i ; $ 2 = <5 2 $2 ■ (2-50) 
Plugging the relations l|2-45l 12-461 12-471 12-50|1 into the expression (|2-43fl for the free energy we obtain 

T=5 4 T + 6 i t 2 \n5 = 5 i F + t 2 ln(5) (2-51) 

Applying the operator Sdg\s=i to both sides we obtain 



= AT- 



J2(2 - K)u K d UK + ^(2 - J)v,,d VJ + 2td t 



K 



T + t 2 , 



(2-52) 



whence the statement of the corollary. Q.E.D. 

We also add a few remarks before moving on with the proof of Thm. 12.11 

Remark 2.2 Formulas (12-421 l^llfll are symmetric in the roles of P and Q: the only non-immediately symmetric term is 
res P 2 QdQ but an integration by parts restores the symmetry. 

°°Q 

Remark 2.3 The formula is derived for polynomial potentials but it could possibly be extended to convergent or formal 
series. 

Remark 2.4 The genus 1 correction to the above formula has been computed in 1141 and is given by 



^ w = -^ln(7 4 £>) 



24 



(2-53) 
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(2-54) 
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Proof of Theorem 12.11 The equivalence of lines (|2-42[) and l|2-43[l follows from the definition of the exterior 
potentials and Lemma |2~T1 

Let us consider the derivative dx '■= g^— done at argument (P or Q) fixed, where the value being kept fixed under 
differentiation is denoted by the corresponding subscript 



Let us set 



(d K $i) Q = f (d K P) Q dQ - P(X q )d K (Q(X q )) 

(d K t> 2 ) P = [ (d K Q)pdP - Q(X p )d K (P(X p )) 
Jx p 

4inF Q := I P$ 1 (Q)dQ + I Q$ 2 (P)dP + - I P 2 QdQ , 

J ooo J oop J con 



' OOQ -/OOp *" J OCQ 

and study the variation of this functional. First off we have the formulas 



oo 



P = V{{Q) - ^ + E KUkQ-x- 1 

OO 

$1 = - Uk ®~ 



K=l 



K=l 



Q = v±(p)-- + Y j JVjP- j - 1 



j=i 



$ 2 = -J^VjP- 



(2-55) 
(2-56) 

(2-57) 

(2-58) 
(2-59) 
(2-60) 
(2-61) 



j=i 



where the Uk , Vj's have been defined in (|2-15|l . We also need the variations of the functions P, Q given here below 

(d K P) Q = Q K - X + 0(Q- 2 ) ; (d K Q) P = 0{P- 2 ) . (2-62) 

With these formulae we can now compute the variation of J-q (the subscript on the loop integrals that follow specify 
the points around which we circulate): 



=.2inU K 



=0 



4ind K F = I (d K P) Q $idQ+ I P {d K $i) Q dQ +<f {d K Q) P $ 2 dP + <f Q (d K $ 2 ) P dP + I P(d K P) Q Q$®^) 

J ooo J ooo Joop J oo p J ooo 



3Q 

= AinUk + <f P f (d k P) Q dQ 
-2nrt [P{X q )d K {Q{X q )) + Q(X p )d K (P(X p ))} 



dQ+ & Q 



(d k Q)pdP 



dP + I P(d K P) Q QdQ 

J OOQ 



AiirUk + <f P 



P 

OOQ 



{d k P) Q dQ 



dQ+ i Q 



(d k Q)pdP 



x„ 



dP + 



[ (d k P) Q dQ 


dQ - I Q 


Ix 


Jx q 


J OOQ 





{d k P) Q dQ 



-2l7Tt 



P(X q )d K (Q(X q )) + Q(X p )d K (P(X p )) - I (d K Q)pdP 



dp 



x„ 



where 



* = 2inU k + 



I p 


l (d k P) Q dQ 


1 OOQ 


Jx q 



dQ-2iirtP{X q )d K {Q{X q )) 



(2-64) 
(2-65) 
(2-66) 

(2-67) 

(2-68) 

(2-69) 



We now use the "thermodynamic identity" 



(d K P) Q dQ = ~(d K Q) P dP , 

so that the last term on line l|2-67|l reads (notice the double change of sign due to the definition of the circle ar 
ooq, which is (homologically) minus the circle around oop) 



OOQ 



f (d k p) Q dQ 


dP = I Q 


/ (d k Q)pdP 


JX q 


J OOp 


Jx q 



dp 



Plugging into the variation of !Fq we get 



P{X q )d K {Q{X q )) + Q{X p )d K (P(X p )) 



(d K Q)pdP 



Finally we claim that the term in the square brackets in lj2-72|) is 



P(X q )d(Q(X q )) + Q(X p )d(P(X p )) - / (dQ)pdP = d 



x„ 



X ~ v 



Q(X P )P(X P 



x„ 



PdQ 



x„ 



PQ-V l {Q)-V 2 {P)) +2<ln( 7 ) 
/ o 

where d denotes any vector field in the space of parameters uk, vj (or even t). Indeed 



Q{X P )P(X P ) 



pX q 




/ PdQ 


= d 


lx p 





X a 



Q(X P )P(X P ) + / P(X)Q'(X)dX 



d(Q(X p ))P(X p ) + Q(X p )d(P(X p )) 



x a 



[(dP) x (X)Q'(X) + P(X)(dQ') x (X)]dX - P{X p )Q'{X p )dX p + P(X q )Q'(X q )dX q = 



X'u 



d(Q(X p ))P(X p ) + Q(X p )d(P(X p )) + 

fX„ 



x„ 



[(dP) x (X)Q'(X) - P'(X)(dQ) x (X)]dX + P{X q )d{Q) x {X q ) - P{X p )d{Q)x{X p ) 

-P(X p )Q'(X p )dX p + P{X q )Q'(X q )dX q . 
In order to proceed we note that if X is a point depending on the parameters, we have 

d(Q(X)) = (dQ)x(X) + Q'(X)dX 

Using this formula we get 



12301 = Q(x p )d(p{x p )) + d{Q(x p ))p{x p ) 

- Q(X p )d(P(X p )) + d{Q(X p ))P(X p ) 



x„ 



x^ p 
x a 



[(dP)x(X)Q'(X) - P'(X)(dQ)x(\)]dX = 



X ,, 



(0P) Q dQ 



which is the term in square brackets in l|2-72l) . Using then Lemma 12. II we have proven that 

2d K T = 2U K - td K {[-Vl(Q) - V 2 (P) + QP] + 2t ln( 7 )} = 2U K - td Kt i . 

The second term is exactly the opposite of the variation of the last two terms in formula (|2-43|l , and this concludes 
the proof (the derivatives w.r.t. vj are treated by interchanging the roles of P and Q). Q.E.D. 

The derivative w.r.t the number operator t has to be treated in a separate way. 
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Corollary 2.2 The derivative of the free energy w.r.t t is given by the formula 

d t T = Q{X P )P{X P ) + [ QdP = 11 . 

Jx p 

Proof. We start with the following observation 



(d t P) Q dQ = -(d t Q)pdP 



dX 
T 



(2-85) 



(2-86) 



In other words they are differential of the third kind with poles at ooq and oop and opposite residues, to wit (using 
the thermodynamical identity l|2-7U|) l 



- idP + 0(P- 2 ) = (d t Q)pdP = -(d t P) Q dQ = idg + 0(Q- 2 ) 
and the differentials have no other singularities. We then proceed as in the proof of Thm. 12.11 





= HQ) 


+ [ (d t P) Q dQ - 
Jx q 


P{X q )d t (Q{X q )) 




= HP) 


+ f (d t Q)pdP- 
Jx p 


Q(X p )d t (P(X p )) = 




= H p ) 
= HP) 


+ f (d t Q) P dP + 
Jx q 

+ f (d t Q)pdP- 
Jx q 


f \d t Q) P dP - Q(X p )d t (P{X p )) 
Jx p 

[ (d t P) Q dQ - Q(X p )d t (P(X p )) 
Jx v 


(d t P) Q 


l 

= ~0 + 


0(Q- 2 ) 




(d t Q)p 


1 

= -p + 


o(p- 2 ) . 





(2-87) 

(2-88) 
(2-89) 
(2-90) 

(2-91) 
(2-92) 
(2-93) 



We then find that 



J OOQ 



=-d\/X 

HQ) + I (WodQ 



dQ + f Q 

OOp 



=d\/X 

HP) + I (d7Q)7dP 



A' 



dP 



—2iirtdt 

= —2iTrtdt 
We now claim that 



Q(X P )P(X P ) 
Q(X P )P(X P ) - 



A ; 
A, 

X„ 



QdP 
QdP 



=dA/A 

PQ (d^PhdQ 



(2-94) 



- 2^7^( J PQ)o + f P [HQ) - ln(A)] dQ + d> Q [ln(P) + ln(A|Ed9B) 

OOQ J OOp 



P [HQ) - ln(A)] dQ = V r 1 (Q) - t Hi) 



and the symmetric formula for the other term. Indeed 



In 



\n-f + 0(Q 1 ) , near ooq 



(2-96) 



(2-97) 
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and thus 



I P [HQ) - ln(A)] dQ= <f (v{{Q) - 1 + 0(Q- 2 )) [hi(Q) - ln(A)] dQ = 

(V{(Q) + 0(Q- 2 )) [HQ) - ln(A)] dQ + 2i7riln( 7 ) 

2ijrt 111(7) = 



(y 1 (Q) + G(Q- 1 )) 

OOQ 

2«7r(F 1 (Q)) + 2 i ^< ln( 7 ) • 



dQ _ dA 
~Q~ ~ A 



Plugging this into l|2-95|) we obtain 



2d t F = -td t 



Q(X P )P(X P ) 



QdP 



X ,, 



Q(X P )P(X P ) 



QdP 



so that 



Q.E.D. 



2dt T = 2 



Q{X P )P{X P ) 



QdP 



(2-98) 

(2-99) 

(2-100) 
(2-101) 

(2-102) 

(2-103) 



2.2 Canonical transformations 

In the Poisson structure l|2-3|) the coordinates ln(A) and t are canonically conjugate, as well as the functions P, Q. 
Therefore it makes sense to find the generating function of these transformations. This was accomplished in the 
context of conformal maps in |.'i0| but it is probably a new result in this context, since now we can express explicitly 
the generating function as integrals of the dToda operators P. Q. I recall that we are looking for a function S!(Q, 
with the property that 

d Qlt fl = PdQ + ln(A)di . (2-104) 
The following proposition gives a representation of such function 

Proposition 2.1 The generating function of the canonical change of coordinates (ln(A),t) 1— > (P, Q) is given by 



fi= / pd Q-Z 

J X Q — 



PdQ-- 



Q{X P )P(X P ) + 
J PdQ + \{Vi{Q) + V 2 (P) -PQ) -t Hi) 



PdQ 



x n 



is the generating function of the canonical transformation (ln(A),t) — > (P, Q), or, 

d Q , t il = d Q n x dQ + {d t n) Q dt = PdQ + ln(A)df . 

Remark 2.5 The function that we get by interchanging the role of P and Q would generate the transformation (— ln(A), i) 
(Q,P) (or the anti canonical one). 



(2-105) 
(2-106) 

(2-107) 



Proof. The first part is obvious 
As for the second we compute 



dntt = P 



by definition. 



(d t n) Q = d t \ Vi(Q) + tHQ) - ^ + o(Q- v 
^HQ)-\dt^ + o(Q- 1 )^ 

-ln(A)+ln(7)-ia tM + (!?(A- 1 ) . 



(2-108) 

(2-109) 

(2-110) 
(2-111) 
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It is well known in the theory of the dToda tau-function that 

d 2 t T = 8^ = 21^) , (2-112) 
and hence the constant term in the above expression vanishes. On the other hand 

(d t fl) Q = f (d t P) Q dQ - P(X q )d t (Q(X q )) = [ ^ - P(X q )d t (Q(X q j) = (2-113) 
Jx q Jx q * 

= ln(A) - ln(A(X g )) - P(X q )d t (Q(X q )) . (2-114) 
By comparison we conclude not only that <9 t f2 = ln(A) but also 

ln(X(X q )) = -P(X q )d t (Q(X q )) . (2-115) 
The proof of the Lemma is complete. Q.E.D. 

3 Planar Free energy for spectral curves of arbitrary genus 

The situation in the case the spectral curve is of higher genus is only slightly different. We will work with the 
following data: a (smooth) curve S 9 of genus g is assigned with two marked points ooq, oop. On the curve we are 
given two functions P and Q which have the following pole structure; 

1. The function Q has a simple pole at ooq and a pole of degree d 2 at oop. 

2. The function P has a simple pole at cop and a pole of degree d\ at ooq. 

From these data it would follow that P, Q satisfy a polynomial relation but we will not need it for our computations. 
By their definition we have 

di+i , . 

P = E u kQ K ^ - n + °^~ 2 ) =: V ^ -q+ °^~ 2 ) ' ncar °°Q t 3 " 1 ) 

K=l ^ ^ 

d 2 + l 

Q = Y J VjP K - 1 -p+0{Q- 2 )=:VJ L {P)-- + 0{Q- 2 ) , near oop. (3-2) 
J=i 

The fact that the coefficient of the power Q^ 1 or P _1 is the same follows immediately from computing the sum of 
the residues of PdQ (or QdP). The formulas for the coefficients uk, vj, t are the same as in the genus zero case, 
viz 

u K = - res PQ~ K dQ , vj = - res QP~ J dP , t = res PdQ = res QdP . (3-3) 

OOQ OOP OOQ OOP 

Note that the requirement that the curve possesses two meromorphic functions with this pole structure imposes 
strong constraints on the moduli of the curve itself. In fact a Ricmann-Roch argument shows that the moduli space 
of these data is of dimension d\ + d 2 + 3 + g; so far our data show only (d\ + 1) + (d 2 + 1) + 1 parameters and 
therefore we add as parameters the following period integrals which, in the matrix-model literature are referred to 
as filling fractions 

£<:= PdQ ' i = 1 '-"S- ( 3 " 4 ) 

Here we have introduced a symplectic basis {a,, bi}i=i... g in the homology of the curve S ff and the choice of the 
a-cycles over the 6-cycles is purely conventional. 

In this extended setting the free energy is defined by the equations 

d UK F g - - res PQ K dQ , d VJ T a - - res QP J dP , (3-5) 

OOQ OOp 

d u T g = l PdQ =:T l . (3-6) 
Jbj 
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Once more we introduce the exterior potentials by the same requirements as in the genus case; 



1 



$ 2 



2l7T 
1 



Q 



*i = ~7HZ 4> ln I 1 - § ) Pd ® = -Vi{Q)+t\n{Q) + I PdQ = OiQ- 1 ) , near oo Q (3-7) 



P 



In 1 - — QdP = -V 2 (P) +t\n(P) + / QdP = ©(P^ 1 ) , near oo 



P 



(3- 



a,. 



In the loop integral formulas the contours wind around the marked points so as to leave the point where the 
potentials are computed inside, i.e., for instance, Q/Q « 1. As in the genus case the two points X p , X q are 
implicitly defined by the requirement C(local parameter). Recall the definition of the "chemical" potential 



/i - Q(X P )P(X P ) + / PdQ 



A' 



(3-9) 



-A,, 



Theorem 3.1 The free energy over the curve S 9 is given by 

1 9 
2T„ = res P$idQ + res Q$ 2 dP + - res P 2 QdQ + tpi + > ej t 



(3-10) 



Proof. It is essentially the same as in the genus zero case. We start from the same expression used in the proof 
there (|2-57[1 and proceed to the variation w.r.t. uk ■ Following the same steps we obtain the expression 



f (d k p) Q dQ 


dQ+ <f Q 


f (d k Q) P dP 


dP + I QPd 


[ (d k P) Q dQ 


Jx q 


J oop 


J X q 


J OOQ 


JXg 



-2int 



P(X q )d K (Q(X q )) + Q(X p )d K (P(X p )) - / (d K Q) P dP 



X c , 



(3-11) 
(3-12) 



where the subscripts in the contour integrals specify the point around which we are circulating. The integration 
by parts of the third term gives 



Aiird K Fo = MirU k + <h P 

J OOQ 



- $ P 

3Q 



= 4i7T[/fc + <b Q 

I OOp 



[ (d k P) Q dQ 


dQ+ <f Q 


[(dhQ)pdP 


Jx q 


J OOp 


Jx q 



[ (d k P) Q dQ 


dQ - I Q 


[ (d k p) Q dQ 


Jx q 


J OOQ 


Jx q 



f (d k Q) P dP 


dP+ <£ Q 


f {d k Q) P dP 


J X q 


J OOQ 


Jx q 



dp + 

dP — 2iirtdKH = 

dP — 2iirtdKh l 



(3-13) 
(3-14) 
(3-15) 



where we have used the thermodynamical identity (|2-7U[) . Due to the genus of the curve the two contours are not 
homologically opposite and the sum of the two integrals gives finally (using the Ricmann bilinear identity which 
we recall in Appendix 



AiirdxTo = AinU k -2^^ / Pd Q <f (^P)gdQ - I PdQ I {d K P)qdQ 

g 

— 4iirU k — 2i7r^^ e$9#-I\ — I\ Qk^i — 2iirtdKH 

9 

= 4iirU k — 2iir ^""^ c^Ok^i — 2intdK^ 



2intd K ^ = (3-16) 
(3-17) 
(3-18) 
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This implies promptly that 

2d K T g = 2d K T a + eid K Ti + td K » = 2U K , (3-19) 

i 

as desired. The check for the variation w.r.t. the filling fractions is the same by noticing that 

(d et P) Q = 0(Q- 2 ) ; (d et Q) P = 0(P- 2 ) , (3-20) 
near the corresponding marked points. Moreover, by definition 

I (d €j P) Q dQ = 6 ij . (3-21) 

From this, following the same formal steps and using the Riemann bilinear identity one easily finds 

g 

23^0 = - Y - td €i iJL , (3-22) 

j=i 

g 

28 u T g = 2d u T Q + td ui i + r, + e A, r j = 2r * ■ (3-23) 

J=l 

This concludes the proof. Q.E.D. 

Remark 3.1 In the two- matrix model setting the moduli of the spectral curve are fixed uniquely in terms of the potentials 
Vi, V2 by the requirement I\ = 0, i = l,...,g (which is a minimum requirement) 

Remark 3.2 The case of the one matrix model is a subcase of this setting where one of the two potentials is quadratic, 
say V2; in this case the spectral curve is hyperelliptic of genus [(di + 2)/2]. 

As for the previous case some special care must be paid for the derivative w.r.t. t, but the result is the same as in 
the genus zero case and it is contained in the next corollary. 

Corollary 3.1 The chemical potential \i is indeed the derivative d t T g . Moreover we have the formula 

/x = res \Vi(Q) - tln(Q / A)l dS - res \V 2 (P) - t ln(PA)l dS - res PQdS + V e< I dS , (3-24) 

OOQ OOp OOQ ^ 

where dS is the normalized differential of the third kind with poles at 00 p t Q and residues ±1 and the function A is 
the following function (defined up to a multiplicative constant) on the universal covering of the curve with a simple 
zero at ooq and a simple pole at oop 

A := exp [ ( dS) . (3-25) 



Proof. The proof is quite more involved than in genus zero and requires some preparation. 

We introduce the normalized differential of the third kind with simple poles at oop g and residues ±1, 

dS = dS OOQOOpl <p dS = , i = 1, . . .g ; res dS = -1 = - res dS . (3-26) 

Ja, °°« °° p 

This provides a coordinate on the (covering of the) curve as follows 

A := e f ds . (3-27) 

Quite clearly the parameter A is defined up to a multiplicative constant and it is multiplicatively multivalued on the 
curve S 9 (around the 6-cycles). This arbitrariness and multivaluedness will not affect our computations because 
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formula l|3-26|) remains unchanged under the rescaling A i— > <5A. Moreover it has no branch-points at oop : Q, where 
instead it has a simple zero and a simple pole (and no other zeroes or singularities). By the definition then 



dS = dln(A) . 

The rest mimics the genus zero case. We introduce the two points A Pj(J such that 



(3-28) 



res [ -Vt(Q) + tln(A) + / PdQ ) dS = 
°°<? V Ja i J 

res [ -V 2 (P) -tln(A) + / QdP ) dS = 
°° p \ Ja p J 



(3-29) 
(3-30) 

From the definition of the points X p ^ q and A Pj9 and following the same steps that were taken in genus we have 

(3-31) 



x„ 



PdQ = res {Vi (Q) - t ln(Q/A)) dS 

A, °°Q 

QdP = - res (V 2 (P) -t\n(P/\)) dS . 



The next relation is different from the genus zero case: 

f a„ 



Indeed we have 
= 



P{A q )Q(A q )+ / QdP = - res PQdS + V a i dS . 



i=0 



(3-32) 



(3-33) 



PdQ-ilnA dS 



QP - Q(A q )P(A q ) - / QdP — tln(A) - / QdP 



= 2i7T res QPdS" + 2wr 



Q(A q )P{A q ) + / QdP 



oo p \ J A 



dS = (3-34) 



QdP - 1 In A ) dS + ^ Ci f dS ,(3-35) 



where we have used the bilinear Riemann identity as well as the fact that § a dS = and the fact that (by definition 
of the point A p ) the term with the under-brace is residue-free at oop 3 . We can now proceed as in Lemma 12. II and 
obtain the formula 



/i = Q(X P )P(X P ) + PdQ = 



-X- D 



= res \Vi (Q) - t \n(Q/ A)l dS - res \V 2 (P) - t ln(PA)l dS - res PQdS + V a I dS . (3-36) 

OOQ OOp OOQ ^ J b 



The formula l|3-36|l is the equivalent in higher genus of the formula in Lemma l2.ll 
Coming back to the proof of the Corollary we note as in Corollary 12 . 21 that 

(d t P) Q dQ = -(d t Q) P dP = dS 

is the normalized Abelian differential of the third kind with poles at the marked points for we have 



= d t €i 



<f (d t P) Q dQ 



(3-37) 



(3-38) 



3 Note that QdP — tdS is a meromorphic Abelian differential without residues at the poles oop q. 
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Computing the variation of J-q we now obtain 



28 t T = -td t 
We now have 



Q(X P )P(X P ) 



QdP 



res PQdS + res In ( j ) PdQ - res ln(PA)QdP 



In ( ^ ) PdQ = res ( K/(Q) - ^ + 0{Q- 2 ) ) In ( ^ ) dQ 



1-2 



Q 



res (V/(Q) + C(Q- 2 )) In -f dQ - 1 res In 



OOQ 



Q 



OOQ \ A 



A 

Q\ dQ 



res (Vi (Q)dS - t ln(Q/A)dS) 



(3-39) 



(3-40) 

(3-41) 
(3-42) 



where, in the second residue, we can replace dQ/Q by dS because ln(Q/A) = 0(1). A similar argument goes for 
the term involving P so that we finally have 

2d t T = -td t fi - res PQdS" + res (Vx(Q) - 1 ln(Q/A)) dS - res {V 2 (P) - iln(PA)) dS . (3-43) 

OOQ OOQ OOp 

Inserting this into the full expression of T g we obtain 

a 

2d t T g = 2d t F + d t (tfi) + tidtFi = (3-44) 

i=l 

= pi - res PQdS + res {VAQ) - tln(Q/A)) dS* - res (F 2 (P) - tln(PA)) dS + V e ; / dS = (3-45) 

OOQ OOQ OOp f— ' J. 

z— 1 1 

= 2m , (3-46) 
where we have used the expression (13- 3611 . This concludes the proof of the corollary. Q.E.D. 

Remark 3.3 The formula for fi seems not symmetric in the roles of P and Q only superficially. In fact, if we exchanged 
the roles, the 3 rd kind differential should also change sign. 

We now investigate the scaling property of this Free energy. First off we have the simple 

Lemma 3.1 Under the change of scale for the functions Q — 8Q and P — oP the chemical potential rescales with an 
anomaly as follows 

H = 8a]l+ Sat ln(Ser) . (3-47) 

Proof. The proof is almost immediate from the expression (|3-24() considering the fact that under that rescaling 
we have 

Uk = <j5 1 ~ k uk ; vj = Sa 1 ~ J vj ; t = Sat , = Sati . (3-48) 
On the other hand the differential dS (and hence the function A) are invariant as follows from its expression 

dS = (d t P) Q dQ = fyPdQ . (3-49) 

The proof follows then immediately from (|3-24[l . Q. E. D. 

With the above lemma we can immediately find the scaling properties of the genus g free energy 

Corollary 3.2 [Scaling properties] The Free energy of the genus g data above satisfies the scaling constraints 

9 



2T g = 



2T a 



L K 



^(1 - K)u K d UK + ^2 v J a vj + ^2 Cid ei + td t 

K J i=l 

9 



K 



T - -t 2 
T a - -t 2 



(3-50) 
(3-51) 
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Proof. The proof is immediate from the definitions of the various objects and using the anomalous scaling of the 
chemical potential fi in Lemma |3. II The constraint i|3-5U|l is obtained by keeping a = 1 and differentiating w.r.t. 
6 at 5 = 1, and l|3-51|) is obtained similarly by interchanging the roles of S and a in the above procedure. Q.E.D. 

The two scaling constraints (|3-5UI 13-5111 form the "two halves" of the scaling 



5^(2 - K)u K d UK + J^i 2 ' J)vjd vj 

. K J 



g 

i=l 



2tdt 



(3-52) 



(obtained by adding 13-501) and l|3-51|l ) which is the translation of the well-known property i|2-44[l for higher genus 
Free energies. Of course the same properties l|3-5(JI l3~3T)l hold also for the free energy in Section |2~T1 fin which case 
the part involving the filling fraction would be missing). 

We conclude with the remark that -quite clearly- the (multivalued) function A is playing essentially the same 
role of the uniformizing parameter in genus zero. The quantities 



ln( 7 ) := - res ln(Q/A)dS , 



ln( 7 ) 



res In(PA)dS , 

OO p 



(3-53) 



are the translation in this setting of the homonymous quantity in the genus zero case, except that they need not 
be equal. However, since now A = exp J dS is defined up to a multiplicative constant depending on the base-point 
of the integral, there would be a choice for the base-point which makes 7 = 7. 



4 Conclusion 

The formulas we have presented fill a gap in both the theory of the dispersionless Toda hierarchy and the two- 
matrix model, where the tau-function (free energy in the planar limit) is known only through its partial derivatives 
but no closed formula for the tau-function itself is known. 

The derivation and the technique of the proof emphasizes the importance of the spectral curve of the model, 
at least in the case of polynomial potentials or, in the dToda language, finite Laurent polynomials for the Lax 
operators. 

On a slightly different perspective, we have computed the free energy of the matrix model in the case where 
the spectral curve is of genus g > 0; the computation is less explicit than in genus zero but the formula is closed. 

It would be interesting to explore further the remnant of the Poisson structure l|2-3|) in this context. In fact it 
is almost immediate to verify that still 

{P, Q} = Xd x Pd t Q - Xd x Qd t P = 1. (4-1) 

The investigation of the Poisson structure will be the topic of subsequent publications. 
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A The Riemann bilinear identity 

It is a classical identity but it may be useful to recall it here. Let be given a curve T, g of genus g and a symplectic 
basis in the homology {a,-, &i}j=i... g . Let 77 and u be two meromorphic Abelian differentials and w be without 
residues and define 

VI = J uj . (A-l) 
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Let S g be a simply connected fundamental domain on the universal covering of the curve then the function f2 is 
single- valued inside E 5 with possibly poles. Let dY> g be the boundary of the domain constituted by the cycles of 
the chosen basis. Then the bilinear Riemann identity claims 



residues 

The proof is very easy and can be found in |12| 
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